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ABSTRACT 

Shafieloo ea al. firstly proposed the possibility that the current cosmic acceleration (CA) is slowing down. 
However, this is rather counterintuitive because a slowing down CA cannot be accommodated in most main¬ 
stream cosmological models. In this work, by exploring the evolutionary trajectories of dark energy equation 
of state wiz) and deceleration parameter qiz), we present a comprehensive investigation on the slowing down 
of CA from both the theoretical and the observational sides. Eor the theoretical side, we study the impact of 
different w(z) by using six parametrization models, and then discuss the effects of spatial curvature. For the 
observational side, we investigate the effects of different type la supernovae (SNe la), different baryon acoustic 
oscillation (BAG), and different cosmic microwave background (CMB) data, respectively. We find that (1) The 
evolution of CA are insensitive to the specific form of w(z); in contrast, a non-flat Universe more favors a slow¬ 
ing down CA than a flat Universe. (2) SNLS3 SNe la datasets favor a slowing down CA at Icr confidence level, 
while JLA SNe la samples prefer an eternal CA; in contrast, the effects of different BAG data are negligible. 
(3) Compared with CMB distance prior data, full CMB data more favor a slowing down CA. (4) Due to the low 
significance, the slowing down of CA is still a theoretical possibility that cannot be confirmed by the current 
observations. 

Subject headings: Cosmology: dark energy, observations, cosmological parameters 


1. INTRGDUCTIGN 


Since its discovery in 1998 ( |Riess et ar]|1998| |Perlmutter| 
|et al.|[l999] ), cosmic acceleration (CA) has been confirmed 
by various cosmological observations, such as type la super- 


novae (SNe la) (Hicken et al. 2009 

Conley et al. 201 ![ Suzuki 

et al. 

2012[ Betoule et al.[[2014[l. 

raryon acoustic oscillation 

(BAG) (Tegmark et al.[[2004[ [Eisenstein et al. 2005 Perci- 

[val et 

al.|2010 [Padmanabhan et al.|2012[), cosmic microwave 

background (CMB) (Spergel et al.|2003[[2007 [Hinshaw et al.[ 

|2013 

Ade et al. |2015a|l, and so 

on. Now CA has become 


)| |Wangl2010[ Li et al.|20iT^ Bamba 

13a| Weinberg et al.|2013| l. Currently, 


& Kamionkowski|2009 

et al.|2012| LI etalW _ 

almost all the mainkream cosmological models, such as A- 
cold-d ark-matter (ACDM) model, holographic dark energy 
model (|^ 2004t [Huang & Li|2004[ Li et al.|2009a b |Zhang| 

|et al.|201^ |Li et al.|2()15b|l, and various scalar fiela dark en- 


ergy ( DE) models Elatev et al.|1999 Armendariz-Picon et al. 


19991 IKamenshchik et al.||2001 

[Caldwell 

12002 1 JPadmanab- 

han|2002 

Wang & Zhang|2008 

Wang et a 

. 2008b, predict an 


etern al CA. _ 

In ( [Shafieloo et al.|[2()09| l, Shafieloo, Sanhi, and Starobin- 
sky firstly proposed the possibility that the current CA is 
slo wing down. By anal ysing the Constitution SNe la sam¬ 
ple (jHicken et al. [200^ with the Chevallier-Polarski-L inder 
(CPU) model ( jChevallier & Polarski|200H|Linder|2003| ), they 
found that the CA has already peaked and is now slowing 
down at \a confidence level (CL). This result is quite in- 
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teresting, because a slowing down CA cannot be accommo¬ 
dated in almost all the m ainstream models, and only some 
strange doomsday models (Wang et al.|2004|l can explain this 


extremely counterintuitive phenomenon. 

In recent years, this topic has attracted a lot of inter- 


ests (|Huang et al. 2009[ |Gong et al.||2010[ |Li, Wu & Yu 

201 1| 

Lin et al.[[20131 

Cardenas &; Rivera||2012l ICardenas 

et al. 

2013 Cardenas 

20141 For example, in (Li, Wu & 

YuP) 

11b, using the C 

^L model, the authors found that the 


ing down CA at Itr CL. In ( [Cardenas & River^|2012| ), also 
using the CPL model, the authors found that the combined 
Union2-(-BAG-(-CMB data also prefer a slowing down CA if 
the spatial curvature is taken into account. In a recent pa¬ 
per, in t he framework of t he CPL model, the Planck Collab¬ 
oration ( [Ade et al.|2015b[ ) reconstructed the evolution of DE 
equation of state (EoS) w{z) by using a combination of Planck, 
galaxy weak lensing and redshift space distortions data; they 
found that wiz) is a decreasing functi on of redshift z when 
z —> 0 (See Fig. 5 of (Ade et al. 2015bll), which also indicates 


a slowing down CA. 

Although the slowing down of CA has been widely inves¬ 
tigated, previous studies mainly focus on the effects of dif¬ 
ferent SNe la datasets; in addition, most of these works only 
consider the CPL model in a fla t Universe. A relative ly com¬ 
prehensive study was given in (Magana et al.[[20T4)l, where 


Magana et al. investigated the evolutionary trajectones of de¬ 
celeration parameter q{z) by using five DE parametrization 
models and four SNe la datasets. In addition, they also stud¬ 
ied the impact of adding CMB and BAG data on the evolution 
of i 7 ( z). However, some other factors are not considered in 
Ref. ([Magana et al.|2014 1. For example, the impacts of spa¬ 


tial c urvature, which may ha ve a degeneracy with the DE EoS 
wiz) ( [Clarkson et al.[|2007| l, are not taken into account. In 
addition, the effects of different BAG and CMB data on the 
evolutionary behavior of CA are not studied, either. To make 
a comprehensive and systematic investigation on this topic, 
all the factors mentioned above need to be taken into account. 
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and this is the main aim of our work. 


For the theoretical side, we study the impact of differ¬ 
ent w{z) by using six popular DE parametrization models, 
including the CPL model, the Jassal-Bagla-Padmanabhan 


(JBP) model i 

J as sal et aL||2005a|b 

1, the Barbosa-Alcaniz 

(BA) model ( 

Barboza & Alcaniz 

20081, the Ma-Zhang 

(MZ) model 

Ma & Zhang [2011 

, the Eeng-Shen-Li-Li 

(FSLL) model 

(jEeng et al. 

20121 and the Wang (WANG) 

model (|Wangj2008 1. In addition, we 

also discuss the effects 


of spatial curvature on this topic. For the observational side, 
we investigate the effects of three kinds of SNe la data, two 
kinds of BAO data, and four kinds of CMB data, respectively. 
The detailed information of these observational data are listed 
in table [T] 

We present our method in Section S our results in Sec¬ 
tion and summarize and conclude in Section 

2. METHODOLOGY 

In this section, we briefly review the theoretical framework 
of various DE models, and describe the observational data 
used in the present work. 


has a dynamical EoS w(z) = wo + Wa^, thus we have 
E(z) = 1 + z)" + 1 +z)^ + O,o( 1 +z)" 

• JBP model Passal et al.|2005a]b| ) has a dynamical EoS 

w{z) = Wo + thus we have 

£(z) = (Oh)( 1 + z)" + O„o( 1 +z)^ + O,o( 1 +zf 

o,, . , \ 1/2 

, ( 6 ) 

• BA model ( |Barboza & Alcaniz|2008| l has a dynamical 

EoS w(z) = Wo + Wa , thus we have 

E(z)= (oH)(l+z)^ + a„o(l+z)^ + 0,o(l+z)2 


2.1. Theoretical Models 

In a non-flat Universe, the Friedmann equation is 


pl^ ~ Pr Pm Pk Pde; ( 1 ) 

where H = d/a is the Hubble parameter, a = (1 is the 
scale factor of the Universe (we take today’s scale factor 
ao= 1), the dot denotes the derivative with respect to cosmic 
time f, Mpi = (SttG)”' is the reduced Planck mass, G is New¬ 
tonian gravitational constant, pr, Pm, Pk and pde are the energy 
densities of radiation, matter, spatial curvature and DE, re¬ 
spectively. The reduced Hubble parameter E{z) = H(z)/Hq 
satisfies 


£2 = U,o(l+z/ + U^o(l+zr + U,o(l+z)" + U^,o/(z), (2) 


where Hq = IQQh km ■ • Mpc~^ is the Hubble constant, 

h is the dimensionless Hubble constant, flro, ^ka and 
UrfeO are the present fractional densities of radiat ion, mat¬ 
ter, spatial c urvature and DE, respectively. Per (|Wang & 
Wang|2013al, we take Uro = Umo/(l +Zeq), where zjq = 2.5 x 
WTuF^„b/2.7K)-4 and = 2.7255K. Since = 
1 - U,„o - U,o - U^o, ^deO is not an independent parameter. 
Here the DE density function /(z) = Pde(z)/Pde(0), which sat¬ 
isfies 

,l+w(zO 


f{z) = exp 


3 / dz' 
Jo 


l+z' 


(3) 


where the EoS w is the ratio of pressure to density for the DE 
W = Pde! Pde- (4) 


MZ model ( Ma & Zhang||20lT] l has a dynamical EoS 
w(z) = wo + wjJ^^^^--Jn2), thus we have 


Eiz) = (u,o(l +z/ + U,„o(l +z)^ + n,o(l+zf 




(8) 


• FSLL model ( |Feng et ar||2012 l has a dynamical EoS 
w(z) = Wo -t Wa jf:p, thus we have 

£(z)= (ao(l+z)'* + U„o(l+z)^ + Uto(l+z)" 

\ 1/2 

• WANG model ( |Wang| |2008| l has a dynamical EoS 
w(z) = wo^ + Wa , thus we have 

£(z)= (Uh)(1+z)" + U„o(1+z)^ + U,o(1+z)2 


For each model, the expression of E{z) will be used to calcu¬ 
late the observational quantities appearing in the next subsec¬ 
tion. To make a comparison, the simplest ACDM model and 
wCDM model are also adopted in the analysis. In addition, 
we also use the deceleration parameter q = to investi¬ 
gate the evolutionary behavior of CA. 


To study the impacts of different w(z), here we consider six 
popular parametrization models: 

• CPL model ( |Chevallier & Polarski|2001| |Linder|2003| l 


2.2. Observational Data 

In this subsection, we describe how to include various ob¬ 
servational data into the statistics. 


“1 In addition to assuming a specific parametrization form for w(z), an¬ 
other popular method of ex ploring the dynamical evolut i on of EoS is adopt- 
ing a binned parametrizatioHHuterer, & Starkman|l20031;|Huterer, & Coorayl 

(20^;|Wang et al.N201U||i()lU;|Li et al.|U(illtV|<iong et al.|U01J|. For' 

simplicity, here we do not study ihe case of binned parametnzation. 


2.2.1. SNe la Data 

Here we use the SNLS3 “Comb ined” samples (consisting 
of 472 SNe la) ( [Conley et al. 2011 1 and the JLA samples (con¬ 
sisting of 740 SNe la) ( jBetoule et al.|2014) l. 
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TABLE 1 

All the observational data used in the present work 


Observation 

Data 

Reference 

SNe la 

SNLS3 (constant ^ 

SNLS3 (linear /3) * 

JLA 

[Conley et al.[[2011[ 

Conley et al. 2011 ; 

Wang & Wang 2013a 

iBetoule et al.|j2014[ 

BAO 

6 dfGS(lD) 

SDSS DR7 (ID) ^ 

BOSSDR9(lD) ^ 

SDSS DR7 (2D) 

BOSS DR9 (2D) 

[Beutler et al.|[201 1[ 

|Padmanabhan et al.|p012J 

|Eisenstein et al.H201 1) 

Hemantha et al. 2014 

Wang 

2014} 

CMB 

Planck 2015 distance priors 

Planck 2013 distance priors 

WMAP9 distance priors 

Planck 2015 full data 

[Ade et al.H2015b[ 

|Wang & Wang|^2013bJ 

|Wang & Wang|^2013bf 

|Aghanim et al.||2015t 


“ “constant /3” means that a constant supernova color-luminosity parameter /3 is used in the analysis. 

* “linear /3” means that a time-varying supernova color-luminosity parameter is used in the analysis. 

“ID” means that the BAO data are obtained by using the spherically averaged one-dimensional (ID) galaxy clustering statistics. 
“2D” means that the BAO data are obtained by using the anisotropic two-dimensional (2D) galaxy clustering statistics. 


The function of th e SNLS3 supernova (SN) data is given 
by ( Conley et al.|201 f) 


XsNLS3 


= C ' AT^, 


( 11 ) 


where AT^ =T^g-T^n,o(i is a vector of model residuals of the 
SN sample, and itib is the rest-frame peak B band magnitude 
of the SN. 

Without considering the possibility of potential SN evolu¬ 
tion, the predicted magnitude of SN can be expressed as 


mmod = 51ogioX>i(z)-a(s- l)-l-^C-l-A^, (12) 


where a and /3 are SN stretch-luminosity parameter and SN 
color-luminosity parameter, s and C are stretch measure and 
color measure for the SN light curve, Af is a parameter rep¬ 
resenting some combination of the absolute magnitude M of 
a fiducial SNe la and the Hubble constant Hq. 

It must be emphasized that, the current studies on var¬ 


ious SNe la data sets, including SNLS3 (fWang & Wang| 

S U nion2.1 (|Mohlabeng & Ralston|[2013| l, and Pan- 
S1 (jScolnic et al.|2014|l, all indicated that, although a 
is still consistent with a constant, should evolve along with 
redshift z at very high confidence level (CL). Moreover, the 
evolution of /3 has significant ef fects on parameter estimation 
of various cosm ological models ( |Wang et al.|2014a|b|c[ 2015[ 
Li et al.||2016||. In addition, it has been proved that ( Wang| 
& Wang|2013a|l the fitting results are insensitive to the func¬ 
tional form of p(z) assumed. So per (Wang et al.|2014a)l, we 


^ It must be emphasized that, in ord er to include host-ga laxy information 
in the cosmological fits, Conley et al. jConley et al.|201l) split the SNLS3 
sample based on host-galaxy stellar mass at 1 O‘'’M 0 , and made At to be dif¬ 
ferent for the two samples. So there are two values of At (i.e. Ati and At 2 ) 
for the SNLS3 data. Moreover, Conley et al. removed Ati and At 2 from 
cosmology-fits by analytically margina lizing over them (for more details, see 
the appendix CofJConleyet^ 201 H ). In the present work, we just follow 
the recipe of ^Conley eTaLj^lTTT , and do not treat At as model parameter. 


also consider the case of a constant a and a linear f3 = /3o+/3iZ, 
then the predicted magnitude of SN becomes 

'«mod = 5 logjo Viiz) - a(s -1) + l3{z)C + M. (13) 

So there are two ways to deal with the SNLS3 data; adopt¬ 
ing a constant /3 = /3o and adopting a linear ^ = /3o + /3iZ. For 
simplicity, hereafter we will call them “SNLS3 (constant /3)” 
and “SNLS3 (linear /3)” data, respectively. 

The luminosity distance I?l(z) is defined as 

VL{z) = Hac-\l+zMz\ (14) 


where c is the speed of light, z and Zhei are the CMB restframe 
and heliocentric redshifts of SN, and the comoving distance 
r(z) is given by 

r(z) = ci/o' |a|-'/2sinn[|H,|'/2r(z)]. (15) 

Here r(z) = sinn(x) = sin(x), x, sinh(x) for fin < 0, 

fin = 0, and fin > 0, respectively. 

The total c ovariance matrix C, which appears in Eq. 11 can 
be written as ([Conley et al.|201 l]l 


C — Dctat “b Cstat “b C 


sys • 


( 16 ) 


Here Dstat denotes the diagonal part of the statistical uncer¬ 
tainty, Cstat and Csys denote the statistical and systematic co- 
variance matrices, respectively. Notice that the covariance 
matrices of the “SNLS3 (constant /3)” and the “SNLS3 (linear 
/?)” data are different. F or the details of const ructing the total 
covariance matrix C, see [Conley et al.|(|201 Ijl. 

In addition, we also use the JLA SN samples ( [Betoule et al. 
2014 1 . Since the function of the JLA samples has the sim- 
ilar form with the case of the “SNLS3 (constant /3)” data, for 
simplicity we do not describe how to include the JLA data 
into the statistics any more. More detail s about the JLA 
samples can be found in ([Betoule et al.[2014|l. 
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2.2.2. BAOData 

Here we use two types of BAO data that are extracted by 
using the spherically averaged ID galaxy clustering statistics 
and the anisotropic 2D galaxy clustering statistics, respec¬ 
tively. For simplicity, hereafter we will call them “BAO(ID)” 
and “BAO(2D)” data, respectively. 

For the BAO observation, two characteristic quantities, 
Dyiz) and rs(zd), are often adopted to constrain various cos¬ 
mological models. The volume-average d effective distance 
Dyiz) is given by (Eisenstein et al. 12005)1, 


Dviz)= 1(1+zrDAz) 


2r, 

HizY ’ 


where the angular diameter distance 

Da(z) = cHqW(z)/(1+z), 


(17) 


(18) 


and the comoving distance r(z) is given by Eq. 15 In addition, 
the comoving sound horizon r^iz) is given by (I Wang & Wang 
12013a| l 


rj(z) = cHq 


-1 


da' 


lo ^/3{1+Rba')a'*EHz'y 


(19) 


where Rb = 31500flboh^(TcMBK)~^, and Elbo is the present 
fractional densit y of baryon. The redshi ft of the drag epoch Zd 
is given in Ref. (|Eisenstein & Hu|1998]l. 

The BAO(ID) data include three data points th at are ex¬ 
tracted from the 6 dFGS ( 6 dF) ( [Beutler et al.|201 1}, the Sloan 


Digital Sky Survey Data Release 7 (SDSS DR7) (Abazajian 
|et al.|2009| l and the Baryon O scillation Spectroscopi c Survey 
Data Release 9 (BOSS DR9) (|Eisenstein et al.|201 l|l, respec¬ 


tively. Mak ing use of the data of 6 df, Beutler et al. (Beutler 
et al!^| 20 lT] l gave 

G(Zrf)/£>v(0.106)^"'“ =0.336 ±0.015. (20) 


Adopting the data of SDSS DR7, Padmanabhan et al. (Pad- 
manabhan et al.| 2012 | l got 

Dv(0.35)/r,(zdY‘“‘‘ = 8.88 ±0.17. ( 21 ) 


Based on th e data of BOSS DR9, Anderson et al. ( [Anderson] 
jet al.| 2012 ] l obtained 


Dv{0.51)/r,{zd)' 


data 


13.67±0.22. 


( 22 ) 


The function of the BAO(ID) data can be expressed as[^ 


In addition, using the 2D t wo point cor relation functions of 
BOSS DR9 samples, Wang f^^ (2014 1 obtained 


//(0.57)r,(zrf)/c‘'“'“ = 0.0444±0.0019, 
D^(0.57)/r,(zrf)‘'“'“ = 9.01±0.23. 

The function of the BAO(2D) data is 

XbAO(2D) = Xl + xi I 
with zi = 0.35 and zi = 0.57. Here 

X? = Aa- [c7'(a,G;)] ^Qh ^Qi = 


(25) 


(26) 


(27) 


where Qx = H( zi)rs{zd)/c, Q 2 = PAiz d/r,( zd), and i = 1 , 2 . 
Based on Refs. ( Hemantha et al.|2014| l and ( |Wang|2014) l, we 
get 

' 0.00000324 -0.00010728 


Ci = 


C 2 : 


-0.00010728 0.0625 


0.00000361 0.0000176111 


0.0000176111 


0.0529 


(28) 


(29) 


2.2.3. CMBData 
For CMB data, we use the distance priors data extracted 


from Planck 2015 (jAghanim et al.|2015 [Ade et al.|2015a|b I, 
Planck 2013 (jWang & Wang||2013bi and WMAP9 (Wang & 


Wang 

the fa 


|2013b|l 

11 ''Plan 


samples. To rnake a comparison, we also use 


Pla nckTT + lowP" data given by the Planck 2015 
data release (jAghanim et al.||2015]l. For simplicity, hereafter 
we will call them “Planck2015”, “Planck2013”, “WMAP9” 
and “Planck2015(full data)”, respectively. 

CMB gives us the comoving distance to the photon¬ 
decoupling surface r(z*) and the comoving sound horizon at 
photon-decoupling epoch rs{z*). The redshift of the photon - 
decoupling epoch z* i s given in Ref. (|Hu & Sugiyama|199^ . 
Wang and Mukherjee ( jWang & Mukherjee|2007| l have shown 
that the CMB shift parameters 


/a = 7rr(z,)/r^(z*), 

R=^Jn,„oH^r(z*)/c, 


(30) 


together with ujb = Etboli^, provide an efficient summary of 
CMB data. Therefore, the x^ function for the CMB distance 
prior data can be expressed as 


3 

Xbao(id) - I 
!=1 




(23) 


Here qi = rs(zrf)/Dv(0.106), q 2 = Dy{0.35)/r^(zd), and ^3 = 
Dv(0.57)/rs(Z(/). In addition, cr, denotes the standard devia¬ 
tion of data, and i= 1,2,3. 

For the BAO(2D) data, we use two observed quantities, 
H{z)rsizd)/c and DA{z)/rs{Zd), that are extracted from SDSS 
DR7 and BOSS DR9 data, respectively. Making use of the 
2D matter power spectr um of SDSS DR7 sam ples, Heman¬ 
tha, Wang, and Chuang ([Hemantha et al.|2014 1 got 


H(0.35)r,(zrf)/c“'“'“ = 0.0431 ± 0.0018, 
DA(0.35)/r,(zrf)"'“'“ = 6.48± 0.25. (24) 


^ Since the redshifts of these three surveys are separated widely, it is rea¬ 
sonable to neglect all the correlations among these surveys. 


xImb = ^Pi [CovcU(P< > Pj)] Ap;, Ap,- = Pi - pf , (31) 

where p\ = R{z*), pi = Liz*), and p^ = ooi,- The covariance 
matrix for (pi,p 2 ,P 3 ) is given by 

CovcMBiPi ,Pj) = cr(pi) a(pj) NoimCovcMBiPi ,Pj), (32) 

where a{pi) is the Icr error of observed quantity /?,, 
NormCowcMBiPi,Pj) is the corresponding normalized covari¬ 
ance matrix, and i,j= 1,2,3. 

As mentioned above, in the present work we use three kinds 
of CMB distance prior data: 

• The Planck 2015 data are ( |Ade et al.|2015b| l 

pf'“ = 1.7382 ±0.0088, 

= 301.63 ±0.15, 

pf'" = 0.02262 ± 0.00029. (33) 
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In addition, 


3.1.1. Impacts of Different w(z) 


NormCovpianckiois = 


1.00 0.64 -0.75 \ 

0.64 1.00 -0.55 

-0.75 -0.55 1.00 / 


(34) 


• The Planck 2013 data are ( |Wang & Wang|2013b[ ) 


pf'“ = 301.57±0.18, 
1.7407 ± 0.0094, 

= 0.02228 ±0.00030. 

In addition, 


NormCoVHa„c,t 2013 = 


1.00 

0.5250 

-0.4235 


0.5250 

1.00 

-0.6925 


-0.4235 

-0.6925 

1.00 


• The WMAP9 data are ( |Wang & Wang|2013b| l 

= 302.02 ±0.66, 

1.7327 ±0.0164, 

= 0.02260 ±0.00053. 

In addition. 


NormCovn^MAP9 = 


1.00 0.3883 -0.6089 \ 

0.3883 1.00 -0.5239 

-0.6089 -0.5239 1.00 J 


(35) 


(36) 


(37) 


(38) 


To make a comparison, we also use the full “Pl anckTT + 
lowP” data given by the Planck 2015 data release (|Aghanim 
et akl |2015|l The complete likelihood description of 


“PlanckiO 15(full data)” can be found in Ref. (|Aghanim et al. 

oft 


2015| l, and the corresponding likelihood softwares can be 
downloaded at http;//pla.esac.esa.int/pla/#results. 


2.2.4. Total function 


The total function is 


- XsN + XbAO + XcMB ■ 


(39) 


In the present work, we mainly make use of the 
“SNLS3(linear ,8)+BAO(lD)+Planck2015” data to study the 
effects of various theoretical factors on CA. In addition, other 
observational data listed in tableware also used for compar¬ 
ison. Finally, we perform an MCMC likelihood anal ysis us¬ 
ing the “CosmoMC” package ( Lewis & Bridle|2002| l. In the 
MCMC analysis, the prior range of each free parameter is uni¬ 
form for all the DE models. These prior ranges are listed in 
Table |2] 


3. RESULTS 

In this section, by exploring the evolutionary trajectories of 
w{z) and q{z), we study the slowing down of CA from both 
the theoretical and the observational sides. 

3.1. Effects of Different Theoretical Models on CA 

Here we discuss the effects of different w{z) and spatial 
curvature on the evolutionary behavior of CA. To study the 
impacts of these theoretical factors, it is necessary to fix 
the observational data used in the analysis first. As men¬ 
tioned above, in this subsection we always make use of the 
“SNLS3(linear ^)-l-BAO(lD)-l-Planck2015” data. 


Eor the theoretical side, let us discuss the impacts of dif¬ 
ferent w(z) first. Table summarizes the fitting results of 
the six DE parametrization models. To make a comparison, 
we also list the fitting results of the ACDM model and the 
wCDM model. Notice that a non-flat Universe is adopted in 
the background. Erom this table we see that, all the six DE 
parametrization models give a similar value of Xmitc which 
is smaller that the Xmin values of the ACDM and the wCDM 
model by ~ 2.5. Besides, for these six models, the differ¬ 
ences of various cosmological parameters (such as U,„o, 
and h) are very small. This means that current observations 
can not distinguish these six dynamical DE models, and thus 
can not put tight constraints on the evolution of EoS w{z). We 
also find that, for the ACDM and wCDM models, the curva¬ 
ture component = 0 (i.e., consistent with a flat Universe) 


at 1(7 CL, which is consistent with the results of many previ- 

ous studies dPadmanabhan et al.|2012 

[Anderson et al. 


2014a|bl [Sanchez et al.||20141 jSamus 

lia et al.[[2014[ [Kazin[ 


However, after taking into ac- 
the curvature component Clko 
will deviate from 0 at 2cr CL. This implies that the spatial 
curvature has a degeneracy with the evolution of EoS w(z), 
and then shows the importance of considering the effects of 
spatial curvature on CA. In addition, we And that, all the six 
dynamical DE models give a current EoS value wq deviating 
from -1 at Icr CL. Moreover, the MZ model give a Wa larger 
than 0 at 1(7 CL, while other five parametrization models give 
a Wa smaller than 0 at 1(7 CL; as will be seen below, these 
constraints on w^ will lead to a slowing down CA. H 

Now we discuss this issue with more details. In Fig.[Tj we 
plot the best-fit results, the Icr and the 2a confidence regions 
of w(z) for the six DE parametrization models. We And that, 
for all the parametrization models, EoS w deviates from a cos¬ 
mological constant at Icr CL, which shows the importance of 
considering the evolution of EoS. Moreover, although differ¬ 
ent DE model gives a different evolutionary behavior of w{z) 
at high redshift, all the six DE models give an increasing w{z) 
when z —0 at almost 2cr CL. This indicates that all the dy¬ 
namical DE models favor a slowing down CA. 

Moreover, we plot the best-fit results, the \a and the 2a 
confidence regions of q{z) for the six DE parametrization 
models in Eig. To make a comparison, we also give the 
corresponding results of the ACDM and the wCDM model. 
It is clear that the ACDM and the wCDM model give a de¬ 
creasing qiz) when z —0, which means that these two models 
prefer an eternal CA. In contrast, for the six DE parametriza¬ 
tion models, q{z) achieves its bottom at redshift z 0.26, and 
then becomes a increasing function of z at more than Icr CL. 
This implies that all the dynamical DE models favor a slowing 
down CA. 

To give a quantitative description on the transition from an 
eternal CA phase to a slowing down CA phase. Here we use 
the transition redshift Zp where the deceleration parameter q{z) 
achieve its bottom. In Table we list the best-fit values of Zp 
given by the six DE parametrization models. It can be seen 
that all the six DE model undergo a phase transition at the 
redshift region [0.23,0.29]; besides, the ESLL model gives a 
largest Zp, while the WANG model gives a smallest Zp. 

Moreover, to assess the statistical significance of the slow- 


’ For the MZ model, is a decreasing function of redshift z, so a Wa 

larger than 0 will lead to a slowing down CA. 
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TABLE 2 

THE PRIOR RANGES OF ALL THE FREE PARAMETERS. THESE PRIOR RANGES ARE UNIFORM FOR ALL THE MODELS CONSIDERED IN THIS WORK. 


Parameter 

a 

/3o 

/3i 


12 60 


h 

Wo 

Wa 

Prior Range 

[0,5] 

[0,5] 

[0,20] 

[-1,1] 

[0,1] 

[0,1] 

[0,5] 

[-6,5] 

[-15,15] “ 


“ For the MZ model, the prior range of Wa is [—35,35]. 


TABLE 3 

EITTING results of THE ACDM, THE wCDM AND THE SIX DE PARAMETRIZATION MODELS. BOTH THE BEST-FIT VALUES, THE Icr ERRORS OF 

VARIOUS PARAMETERS AND THE xL,i OF VARIOUS MODELS ARE LISTED. 


Parameter 

ACDM 

wCDM 

CPL 

JBP 

BA 

MZ 

FSLL 

WANG 

a. 

fo 

01 

12*0 

1260 

12 m 0 

h 

Wo 

Wa 

Sp- ■ 

■^min 

1 4^4+0.098 

I-^^o_q.360 

-c o 1 q+0.978 

-0.0009" 

0-0483)|«5f 

0.300i«l« 

0.684f»;[]]« 

386.5 

1 4'2144).098 

111 

1 469+^'^’^^ 

1 

O.U/U_i 126 

-0.0013" 

0.0480^;“P 

0.297t«[| 

0.687" 

-1 ni4+^'^^o 
i.UiH-ooso 

386.5 

1 400 + 0.100 

4-4AD_0,374 
c oqc+1.009 

152 

U.ZVO_0015 

0 

-0705» 

-2.286A;fJ 

383.7 

1 qq-;v-K)-099 

4-a:?j_0_io8 

4-a^a_0,368 

4.668fl»| 

-0 0099+^4039 

oZo^f 

0.295«;||l^| 

0.687" 

-0.648):*4)2 

-3 4191^:^^^ 

286 

384.1 

1 72Q7+0.1UJ 
1.35/_0 103 

1-437" 

5.085fl:Op 

U.U‘t/A_onQ24 

0 994+0-^13 
u.zy^_00i4 

ooiVm 

U.80 0,176 
-1 

i.z,uu_og 8 o 

383.8 

1 400+0.098 

i.‘fUy_o 110 

1 QQ4+0.410 
i-A^‘t_0.368 

C TTO+O.^OS 

153 

'-'•^I4 1_oqq52 
U.U^t/OnoQJS 
0.295tg» 
0.689" 

-0 764+®-'“ 
«405)ff 
383.8 

1 '1Q7+U.U9/ 

^•^^'-0.109 

1 44q+0,409 

5.035!?:?«| 

-0-0118S 

u.u^oy_0.QQ24 

0 900+0.013 

U.ZVU_ooi3 

U.D^A_0016 
-0 690^^-^^^ 

Kj.{jyv_Q2QA 

7 t'jc+1.291 
z,.jjj_i_27g 

383.7 

1 447+0.099 

1 77 C+ 0.406 
^•■^'‘^-0.371 
c 107+1.008 
^•^^^-1,134 

0-0482)|“| 

0-299i“l! 

0.685" 

-0 

-Q-2Q5 
-1 446+^-3^^ 

4 •‘+^o_0.292 

384.0 


TABLE 4 

The Xmin®’ ^AICS AND ABICS OF THE ACDM AND THE SIX PARAMETRIZATION MODELS. 


Parameter 

ACDM 

CPL 

JBP 

BA 

MZ 

FSLL 

WANG 

^miii 

386.5 

383.7 

384.1 

383.3 

383.8 

383.7 

384.0 

AAIC 

0 

1.2 

1.6 

1.3 

1.3 

1.2 

1.5 

ABIC 

0 

9.5 

9.9 

9.6 

9.6 

9.5 

9.8 


TABLE 5 

Best-fit results of transition redshift Zp given by the six parametrization models. 


Pai'ameter 

CPL 

JBP 

BA 

MZ 

FSLL 

WANG 

Zp 

0.258 

0.241 

0.269 

0.254 

0.286 

0.238 


ing down tr end of CA, we adopt the Akaike information cri¬ 
teria ( AIC) (|Akaike| \ 91A\ and Bayesian information criteria 
(BIC) ( |Sch\varz|1978| l that are defined as 


AIC = Xmin + 2^, BIC = Xmin + ^ In 


(40) 


where k is the number of free parameters, and N is the num¬ 
ber of total data points. In Table we list the results of 
Xm,„, AAICmodei and ABICmodei for the ACDM and the six 
DE parametrization models, where 


A AlCmodel — AlCmodel AIC ACDM; 

ABICtnodel = BlCmodel-BIC ACDM- (41) 


From this table we see that, compared with the ACDM model, 
all the six DE parametrization models yield larger AICs and 
BICs. This means that the dynamical evolution of EoS is not 
favored by the current observational data. 

In conclusion, once the dynamical evolution of DE EoS 
w is taken into account, the combined SNLS3(linear /3) + 
BAO(1D)-I-Planck2015 data prefer a slowing down CA at 
more than la CL. Since all the six DE models give same 
evolutionary trends of w(z) and q(z) when z —> 0, our con¬ 
clusion is insensitive to the specific functional form of w(z). 
This conclusion is consistent with the results of Ref. ( |Magana 
|et al.|[MT^. In addition, due to low statistical signmcance. 


the slowing down of CA is still a theoretical possibility that 
cannot confirmed by the current observations. 

3.1.2. Impacts of Spatial Curvature 

Then, let us turn to the impacts of spatial curvature. Table|g 
shows the fitting results of the CPE, JBP and FSLL models, 
where both the cases with and without spatial curvature are 
taken into account. We find that, in a flat Universe, these 
three parametrization models give a wq consistent with -1 and 
a Wa consistent with 0 at Icr CL. In contrast, in a non-flat Uni¬ 
verse, Wo deviates from -1 at Icr CL, while Wa deviates from 
0 at Icr CL. This means that there is a significant degener¬ 
acy between the spatial curvature and the EoS w(z), which is 


consistent with the results of ([Clarkson et al. 20071 

Anderson] 

|et al.|2014a|bl|Sanchez et al.|2014[|Kazin et a .|2014 

L Inaddi- 


tion. It IS clear that this conclusion is insensitive to the specific 
functional form of w(z). These results show the importance of 
considering spatial curvature in the study of CA. 

In Fig. 1^ we plot the Icr and 2cr probability contours in 
the wo-Wfl plane (left panel) and the Icr confidence regions of 
w(z) (right panel) for the CPE model. From the left panel we 
see that, the CPE model is consistent with the ACDM model 
in a flat Universe. In contrast, the CPE model deviates from 
the ACDM model at Icr CL in a non-flat Universe; in other 
words, in the framework of non-flat Universe, the CPE model 
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Fig. 1. — (color online). The best-fit results (red dash-dotted line), the \<j confidence regions (black solid lines) and the 2c7 confidence regions (green dashed 
lines) of EoS w(z) for the six DE pai'ametrization models. 


TABLE 6 

Fitting results of the CPL, the JBP and the FSLL models, where both the best-fit values and the la errors of various parameters 

ARE LISTED. THE LEGENDS “ WITH 0/.o” AND “ WITHOUT REPRESENT THE CASES WITH AND WITHOUT CONSIDERING SPATIAL CURVATURE, 

respectively. 




CPL 


JBP 

FSLL 

Parameter 

WITH \ li(Q 

WITHOU I ^IkQ 

WITH i 

WITHOUT Ojto 

WITH 12*0 

WITHOUT S 2*0 

a 

A 

/9i 

^hO 

^mO 

h 

Wo 

Wa 

■ 

1 400 + 0.100 

i.‘fOV_Q 1 Q 9 

i.jjo_Q_374 

C ^QC + l.OO? 

-o.oiostJSi^ 

0.295j|p 

0.689;i|:l 

-0.705^:™ 

-2-286fl:“ 

383.7 

1 421+0.099 

1.4tOO_Q 37 ^ 
c 19 Q+I.OOI 

u.U4/y n.p24 

0,300!j;g 

0 , 686 ti 

-0 QS2+0134 

u.uoz ,, 441, 

386.6 

1 TjQc+o.oyy 

4.668!|“f 

-0.0099!«;“9 

0-0480!«;«<>P 

0.295!flJ[| 

0-6873;i|j5 

-0.648!«|| 

-3 419+2.290 

J-^>-y_2286 

384.1 

L384_q jiQ 

1 4T7+0-T03 

i.‘+j/_0_374 

S 226+^-^^* 

O.ZZD_i 139 

0 0477+^-^^23 

0 298+^-^^^ 

aeLlil 

-0 960+° '*’ 

’’■++”’-0.179 
-0 717+1.341 
-1.149 

386.6 

1 OQ 7 -I-U.O 97 

1 4 40+0-409 

1 .‘+H -^_0 37o 
O.UJO_i 164 
n ni 1 q+0.0038 

U.U118_DQQ5g 

0 0469+^-^^22 

U.U40V_ppp24 

0 290+^-^^^ 

U.Z!^U_Qp ,3 

0.695tl|j5 

-0.690^1'h 

-2-335+1 

383.7 

1 41^+0,099 

1 466+^-^*l 

i.4fDO_o 374 

4 Q7«:+L007 

0 n474+^'®^23 
U.U^/^_ppp25 

0 295+^'^*^ 
069h 8:8I? 
-0.968!°:1« 
-0.165!°;g‘ 
386.6 
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ACDM 







0,4 
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0.2 
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•ii -0.2 



D" 
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-0.4 
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-0.6 


— Best-fit 
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— Iff CL 


-0.8 


--- 2ffCL 


-^■8 




WCDM 





Fig. 2.— (color online). The best-fit results (red dash-dotted line), the Icr confidence regions (black solid lines) and the 2cr confidence regions (green dashed 
lines) of deceleration parameter q{z) for the ACDM model, the wCDM model and the six DE parametrization models. 
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Fig. 3.— (color online). The Icr and 2cr probability contours in the WQ-Wa plane (left panel) and the \<j confidence regions of w(z) (right panel) for the CPL 
model. Black solid lines denote the results obtained in a non-flat Universe, while red dashed lines denote the results obtained in a flat Universe. The legends “ 
WITH n^-o” and “ WITHOUT represent the cases with and without spatial curvature, respectively. In the left panel, to make a comparison, the fixed point 
{wQ^Wa} = {—1,0} of the ACDM model is marked as a cyan star. In addition, the magenta dashed lines divides the panel into two regions: big rip region and no 
big rip region. 


will predict a Universe without big rip. From the right panel 
we find that, in a flat Universe, w(z) is consistent with -1, 
which corresponds to an eternal CA. In contrast, in a non¬ 
flat Universe, w{z) is a decreasing function of redshift z when 
z —0, which corresponds to a slowing down CA. 

Fig. a shows the best-fit results (upper left panel), the Icr 
confidence regions (upper right panel), and the 2a confidence 
regions (lower panel) of deceleration parameter qiz) for the 
CPU model. We see that, in a flat Universe, q{z) is always an 
increasing function of z, which corresponds to an eternal CA 
In a non-flat Universe, qiz) becomes an decreasing function 
when z 0 at Icr CL, but is still consistent with an increasing 
function at 2a CL. This means that the CPL model in a non¬ 
flat Universe favors a slowing down CA at Ict CL. 

In conclusion, considering spatial curvature or not will sig¬ 
nificantly change the evolutionary trajectories of CA. Using 
the SNLS3(linear /3)-l-BAO(lD)-l-Planck2015 data, a non-flat 
Universe prefers a slowing down CA at Icr CL, while a flat 
Universe favors an eternal CA. 

3.2. Ejfects of Different observational Data on CA 

Now, we discuss the impacts of different observations, in¬ 
cluding SNe la, BAO, and CMB data, on the evolutionary be¬ 
havior of CA. To study the impacts of different observational 
data, it is necessary to fix the cosmological model in the back¬ 
ground first. So in this subsection, we always use the CPL 
model in a non-flat Universe. 


suits of the “SNLS3(linear /3)” data. 

In Fig. we plot the best-fit results (upper left panel), the 
IfT confidence regions (upper right panel) and the 2tT confi¬ 
dence regions (lower panel) of deceleration parameter qiz) 
given by three kinds of SNe la data. We see that, both the 
best-fit results of qiz) given by the “JLA”, the “SNLS3(linear 
/3)” and the “SNLS3(constant /?)” data correspond to a slow¬ 
ing down CA. At Icr CL, the results given by the “JLA” data 
are consistent with an eternal CA, while the “SNLS3(linear 
/3)” and the “SNLS3(constant /?)” data still favor a slowing 
down CA. At 2cr CL, both the “JLA”, the “SNLS3(linear /3)” 
and the “SNLS3(constant /3)” data are consistent with an eter¬ 
nal CA. This means that the SNLS3 data favor a slowing down 
CA at Icr CL, while the “JLA” data prefer an an eternal CA at 
Icr CL. 

In conclusion, the use of SNe la data has significant impacts 
on the evolutionary behavior of CA: the SNLS3 datasets favor 
a slowing down CA at Icr CL, while the JLA samples prefer 
an eternal CA. Without considering s patial curvature, Ma gana 
et al. found that the Union2.1 data ( [Suzuki et al.| 2012)1 pre¬ 
fer an eternal CA; i n contrast, the Constituti on ( jFticken et al.j 
|2009| l, the Union2 (jAmanullah et al.||2010 l and the LOSS- 
Union ( jGaneshalingam et al.|2013 1 SNe la d ata favor a slow¬ 
ing down CA at Icr CL ( [Magana et al.|2014 1 . There fore, our 
study verifies the conclusion of[Magana et al.1([2014[l. 


3.2.2. Impacts of Different BAO Data 


3.2.1. Impacts of Different SNe la Data 

For the observational side, let us discuss the impacts of dif¬ 
ferent SNe la data first. In Fig. 1^ we plot the Icr and 2a 
probability contours in the wo-Wa plane (left panel) and the Icr 
confidence regions of w(z) (right panel) given by three kinds 
of SNe la data. From the left panel we see that, the Icr con¬ 
tour given by the “JLA” data accommodates the fixed point of 
the ACDM model; in contrast, the Icr contours given by the 
“SNLS3(linear /3)” and the “SNLS3(constant /?)” data devi¬ 
ate from the result of the ACDM model. From the right panel 
we And that, the Icr region of w(z) given by the “JLA” data 
is consistent with -1, which corresponds to an eternal CA. In 
contrast, the Icr regions of w(z) given by the “SNLS3(linear 
/3)” and the “SNLS3(constant /?)” data are decreasing func¬ 
tions of z when z —0, which corresponds to a slowing down 
CA. In addition, the results given by the “SNLS3(constant /?)” 
data are a little closer to a cosmological constant than the re¬ 


Then, let us turn to the impacts of different BAO data. In 
Fig.l^ we plot the Icr and 2a probability contours in the wq- 
Wa plane (left panel) and the Icr confidence regions of wiz) 
(right panel) given by two kinds of BAO data. From the left 
panel we see that, the Icr contours given by the “BAO(ID)” 
and the “BAO(2D)” data cannot accommodate the fixed point 
of the ACDM model. From the right panel we And that, the 1 cr 
regions of wiz) given by the “BAO(ID)” and the “BAO(2D)” 
data are decreasing functions of z when z —0, which corre¬ 
sponds to a slowing down CA. In addition, the “BAO(ID)” 
data give a little larger qiz), which is a little closer to the case 
of the ACDM model. 

In Fig. we plot the best-fit results (upper left panel), the 
Icr confidence regions (upper right panel) and the 2cr confi¬ 
dence regions (lower panel) of deceleration parameter qiz) 
given by two kinds of BAO data. We see that, for both the 
cases of “BAO(ID)” and “BAO(2D)”, qiz) undergoes a transi¬ 
tion at z — 0.25, an then becomes an decreasing function when 
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Fig. 4.— (color online). The best-fit results (upper left panel), the Itr confidence regions (upper right panel), and the 2a confidence regions (lower panel) of 
deceleration parameter q{z) for the CPL model. Black solid lines denote the results obtained in a non-fiat Universe, while red dashed lines denote the results 
obtained in a fiat Universe. The legends “ WITH and “ WITHOUT Qho” represent the cases with and without spatial curvature, respectively. 




Fig. 5.— (color online). The Icr and 2a probability contours in the w^-Wa plane (left panel) and the Itr confidence regions of w(z) (right panel) given by three 
kinds of SNe la data. The legends “SNLS3(linear /?)”, “SNLS3(constant /3)” and “JLA” represent “SNLS3(linear /3)-l-BAO(lD)-l-Planck2015”, “SNLS3(constant 
/3)-l-BAO(lD)-tPlanck2015” and “JLA-HBAO(lD)-tPlanck2015” data, respectively. Black solid lines denote the results given by the “SNLS3(linear /3)” data, red 
dashed lines denote the results given by the “SNLS3(constant /3)” data, while green dash-dotted lines denote the results given by the “JLA” data. In the left panel, 
the fixed point {wo,Wa} = {—1,0} of the ACDM model is marked as a cyan star. In addition, the magenta dashed lines divides the panel into two regions: big rip 
region and no big rip region. 


z —I 0 at Icr CL; In addition, this feature disappears at 2a CL, 
This shows that both the “BAO(ID)” and the “BAO(2D)” data 
favor a slowing down CA at Icr CL. Moreover, we find that 
the differences between the evolutionary trajectories of q(z) 
given by the “BAO(ID)” and the “BAO(2D)” data are very 
small. 

In conclusion, the effects of different BAO data on CA are 
negligible. 

3.2.3. Impacts of Different CMB Data 

Finally, let us discuss the impacts of different CMB data. 
To do this, we make two comparisons for the CMB data: 1. 
a comparison among different CMB distance prior data; 2. 
a comparison between the CMB distance prior data and the 
full CMB data. It must be emphasized that, in order to get 


better visual effects about these comparisons, in this subsub¬ 
section we make use of the JLA SNe la samples, instead of 
the “SNLS3(linear /3) data. 

In Fig. we plot the Icr and 2cr probability contours in 
the WQ-Wa plane given by three kinds of CMB distance prior 
data (upper left panel) and by two kinds of Planck2015 data 
(upper right panel), as well as the Itr confidence regions of 
w{z) given by three kinds of CMB distance prior data (lower 
left panel) and by two kinds of Planck2015 data (lower right 
panel). As mentioned above, to get better visual effects about 
the comparisons of various CMB data, the JLA data are used 
in the analysis, and this is why the Icr contours given by the 
“Planck2015”, the “Planck2013”, and the “WMAP9” data can 
accommodate the fixed point of the ACDM model (see Fig. 
for details). From the left two panels we see that, both the 
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Fig. 6. — (color online). The best-fit results (upper left panel), the Icr confidence regions (upper right panel) and the 2g confidence regions (lower panel) 
of deceleration parameter q{j) given by three kinds of SNe la data. The legend “SNLS3(linear ^)”, “SNLS3(constant /3)” and “JLA” represent “SNLS3(linear 
/3)-l-BAO(lD)-l-Planck2015”, “SNLS3(constant ^)+BAO(lD)+Planck2015” and “JLA+BAO(lD)-l-Planck2015” data, respectively. Black solid lines denote the 
results given by the “SNLS3(linear /3)” data, red dashed lines denote the results given by the “SNLS3(constant /3)” data, while green dash-dotted lines denote the 
results given by the “JLA” data. 




Fig. 7.— (color online). The Icr and 2cr probability contours in the wq-Wq plane (left panel) and the Icr confidence regions of w{z) (right panel) 
given by two kinds of BAO data. The legend “BAO(ID)” and “BAO(2D)” represent the “SNLS3(linear /3)-l-BAO(lD)-l-Planck2015” and “SNLS3(linear 
^)-i-BAO(2D)-i-Planck2015” data, respectively. Black solid lines denote the results given by the “BAO(ID)” data, while red dashed lines denote the results 
given by the “BAO(2D)” data. In the left panel, the fixed point {wq, Wo} = {—1,0} of the ACDM model is marked as a cyan star. In addition, the magenta dashed 
lines divides the panel into two regions: big rip region and no big rip region. 


probability contours of {wo,Wa} and the confidence regions 
of w{z) given by the three kinds of CMB distance prior data 
are almost overlap; this means that the impacts of different 
CMB distance prior data on CA are negligible. In contrast, 
from the right two panels we hnd that, the Icr contours given 
by the “Planck2015(full data)” data cannot accommodate the 
hxed point of the ACDM model; correspondingly, the Icr re¬ 
gion of w{z) given by the “Planck2015(full data)” data is a 
decreasing function of z when z —> 0, which corresponds to 
a slowing down CA. In addition, the differences between the 
“Planck2015(full data)” results and the “Planck2015” results 
are quite large. 

In Fig. we plot the best-fit results (upper left panel), 
the Ict coimdence regions (upper right panel) and the 2a con¬ 


fidence regions (lower panel) of deceleration parameter q{z) 
given ^ the “Planck2015” and the “Planck2015(full data)” 
data. Notice that the JLA data are used in the analysis 
for better visual effects. From this figure we see that, both 
the best-fit results of ^(z) given by the “Planck2015” and the 
“Planck2015(full data)” data correspond to a slowing down 
CA. At Icr CL, the results given by the “Planck2015” data 
are consistent with an eternal CA, while the “Planck2015(full 
data)” data still favor a slowing down CA. At 2cr CL, both the 
“Planck2015” and the “Planck2015(full data)” data are con- 

* Since the differences among the evolutionary trajectories of q(z) given by 
three kinds of CMB distance prior data are very very small, the corresponding 
results are not discussed here. 
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Fig. 8.— (color online). The best-fit results (upper left panel), the Icr confidence regions (upper right panel) and the 2cr confidence regions (lower panel) of 
deceleration parameter q{z) given by two kinds of BAO data. The legend “BAO(ID)” and “BAO(2D)” represent the “SNLS3(linear /3)+BAO(lD)+Planck2015” 
and “SNLS3(linear /3)-i-BAO(2D)-i-Planck2015” data, respectively. Black solid lines denote the results given by the “BAO(ID)” data, while red dashed lines 
denote the results given by the “BAO(2D)” data. 
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Fig. 9.— (color online). The Icr and 2cr probability contours in the wq-Wo plane given by three kinds of CMB distance prior data (upper left panel) and 
two kinds of Planck2015 data (upper right panel), as well as the Icj confidence regions of w(z) given by three kinds of CMB distance prior data (lower left 
panel) and two kinds of Planck2015 data (lower right panel). The legends “Planck2015”, “Planck2013”, “WMAP9”, and “Planck2015(full data)” represent the 
“JLA+BAO(1D)+Planck2015”, “JLA+BAO(lD)-HPlanck2013”, “JLA-hBAO(1D)-i-WMAP 9”, and “JLA+BAO(1D)+Planck2015(full data)” data, respectively. In 
the left two panels, black solid lines denote the results given by the “Planck2015” data, red dashed lines denote the results given by the “Planck2013” data, while 
green dash-dotted lines denote the results given by the “WMAP9” data. In the right two panels, black solid lines denote the results given by the “Planck2015” 
data, while red dashed lines denote the results given by the “Planck2015(full data)” data. 
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Fig. 10.— (color online). The best-fit results (upper left panel), the Icr confidence regions (upper right panel) and the 2cr confidence regions (lower panel) of 
deceleration parameter q{z) given by the “Planck2015” and the “Planck2015(full data)” data. The legends “Planck2015” and “Planck2015(full data)” represent 
the “JLA-l-BAO(lD)-l-Planck2015” and “JLA+BAO(lD)-l-Planck2015(full data)” data, respectively. Black solid lines denote the results given by the “Planck2015” 
data, while red dashed lines denote the results given by the “Planck2015(full data)” data. 
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sistent with an eternal CA. This means that after adopting the 
JLA data, the “Planck2015” data favor an an eternal CA at Icr 
CL, while the “Planck2015(full data)” data prefer a slowing 
down CA at Icr CL. 

In conclusion, compared with CMB distance prior data, full 
CMB data more favor a slowing down CA. As far as we know, 
the effects of different BAO and CMB data on CA have not 
been studied in the past. 


a theoretical possibility that cannot confirmed by the 
current observations. 

• (3) Considering spatial curvature or not will signifi¬ 
cantly change the evolutionary trajectories of CA: in the 
framework of dynamical DE models, a flat Universe fa¬ 
vors an an eternal CA, while a non-flat Universe prefers 
a slowing down CA (see Table]^ Fig.|^and Fig.Q. 


4. SUMMARY AND DISCUSSION 

As is well known, in recent 15 years, CA has become 
one of the biggest puzzles in modem cosmology. Cur¬ 
rently, almost all the main stream cosmological m odels pre¬ 
dict an eternal CA. But in ( |Shafleloo et al.| 2009| l, Shafleloo, 


Sanhi, and Starobinsky proposed the possibility that the cur¬ 
rent CA is slowing down. Since this extremely counterin¬ 
tuitive phenomenon cannot be accommodated in almost all 


terests (|Huang et al.||2009| |Gong et al.||2010 |Li, Wu & Yu 

2011 

|Lin et al.||2013l ICardenas & Rivera||2012l ICardenas 

et al. 

2013 

Cardenas |2014|l in recent years. 


different SNe la datasets; in addition, most of these works 
only consider the CPU model in a fla t Universe. A rela - 
tively comprehensive study was given in (|Magana et al.|2014|l, 
where Magana et al. investigated the evolutionary trajectories 
of deceleration parameter q{z) by using five DE parametriza- 
tion models and four SNe la datasets. However, some other 
factors, including the impacts of spatial curvature, as well as 
the effe cts of different BAO and CMB data, are not considered 
in Ref. (|Magana et al.|2014|l. 

In the present work, by taking into account all the factors 
mentioned above, we perform a comprehensive and system¬ 
atic investigation on the slowing down of CA from both the 
theoretical and the observational sides. For the theoretical 
side, we study the impact of different w{z) by using six DE 
parametrization models, including the CPL, the IBP, the BA, 
the MZ, the FSLL and the WANG model. In addition, we also 
discuss the effects of spatial curvature on this topic. For the 
observational side, we investigate the effects of three kinds of 
SNe la data, two kinds of BAO data, and four kinds of CMB 
data, respectively. The detailed information of these observa¬ 
tional data are listed in table [1] 

Our conclusions are as follows: 


• (1) The evolutionary behavior of CA are insensitive 
to the specific form of w{z) (see Table Fig. 0 and 
This result i s consistent with the conclusion of 


Magana et al. (2014 1 . 


• (2) Compared with the ACDM model, the dynamical 
evolution of EoS is not favored by the current observa¬ 
tional data (see table |^. This means that, due to low 
statistical significance, the slowing down of CA is still 


(4) The use of SNe la data has significant impacts on 
CA: for the CPL model, SNLS3 datasets favor a slow¬ 
ing down CA at 1(7 CL, while JLA samples prefer an 
eternal CA (see Fig. [5]and Fig.lhll. The se results verify 
the conclusion of Magana et al.]([2014|). 


• (5) The effects of different BAO data on the evolu¬ 
tionary behavior of CA are negligible (see Fig. 1^ and 
Fig.[§. 


• (6) Compared with CMB distance prior data, full CMB 
data more favor a slowing down CA in the framework 
of the CPL model (see Fig|^and Fig. [TOj i. As far as we 
know, the effects of different CMB data on CA have not 
been studied in the past. 


It is clear that the evolutionary behavior of CA depends on 
both the theoretical models and the observational data: for the 
theoretical side, it will be significantly changed by the intro¬ 
duction of spatial curvature; for the observational side, it is 
very sensitive to the SN la and the CMB data. 

Since SNe la data has significant impacts on CA, it is very 
important to take into account the systematic error of SNe 
la data seriously. In the present work, only the evolution of 
SN color-luminosity parameter /3 is considere d. Other factors, 
such as the evolution of intrinsic scatter (Tint (|Marriner et al. 


201 l|l, the different choice of SN lightcurve fitter model (Ben 


gochea & Rossi|[2014 IHu et al.]|2015a|bl l, and the different 
method of calculating SN lightcurve parameters (|Dai & Wang| 


2015|l, can also affect the cosmology-fit results given by SNe 
la data. These topics deserve further investigations in future. 

One of the keys of exploring CA is to break the degeneracy 
between the spatial curvature and the EoS parameters with 
more data. It will be interesting to study the impacts of other 
cosmological observations (such as the Hubble parameter, the 
weak lensing, and the growth factor data) on CA. This issue 
will be studied in a future work. 
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